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Abstract A methodology based on point-like sources is discussed, enabling a reliable 
estimate of the acceleration of the Pioneer 10 and 11 probes caused by thermal effects. 
A sensitivity analysis of the several parameters of the model allows for a clear indication 
of the possible thermal origin of the so-called Pioneer anomaly. 
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1 Introduction 

1.1 General Background 

A decade ago, an anomalous, sunbound acceleration on the Pioneer 10 and 11 probes 
was discovered, using two independent code analyses [1,2]; this acceleration may be 
characterized as constant, with a magnitude of ap; ~ (8.5 ± 1.3) x 10~ 10 m/s 2 . Early 



Dcpartamcnto de Fi'sica, Institute- Superior Tecnico; 
also at Instituto de Plasmas e Fusao Nuclear 

O. Bertolami 

E-mail: orfeu@cosmos.ist.utl.pt 

J. Paramos (Speaker) 

E-mail: jorge.paramos@ist.utl.pt 

Departamento de Engenharia Mecanica, Instituto Superior Tecnico; 
also at Centro de Ciencias e Tecnologias Aeronauticas e Espaciais 

F. Francisco 

E-mail: frederico.francisco@ist.utl.pt 
P. S. J. Gil 

E-mail: p.gil@dem.ist.utl.pt 



Address: Av. Rovisco Pais 1, 1049-001 Lisboa, Portugal 



2 



discussions [1] argued that the effect cannot be explained in terms of a misestimation 
of the systematic effects (of thermal nature, or due to electric or magnetic forces, solar 
radiation and solar wind pressure, mechanical defects or errors in the Doppler tracking 
algorithms used, etc.), despite claims otherwise [3]. 

Amongst these "conventional" candidates, one of the prime contenders is the reac- 
tion force due to thermal radiation from the main bus compartment and the radiother- 
mal generators (RTGs) , either directly pointing away from the Sun, or reflected by the 
main antenna dish. Quite obviously, an acceleration due to thermal dissipation should 
present a secular evolution parallel to the available power in the RTGs, thus decaying 
at an approximately exponential rate. This said, one must note that Ref. [2] shows that 
the available data may be fitted to such a signature, i.e. a linear decay with a time 
constant larger than 70 years: since the half-life of the plutonium source in the RTGs 
is ~ 88 years (that should be effectively less due to thermal coupling degradation) , this 
indicates that thermal effects could account for the Pioneer anomaly. This scenario 
is being extensively scrutinized by groups within the Pioneer collaboration team [4, 
5], and has been the motivation of the recent study carried out by our group [6]. A 
consistent study of secular and spatial trends is of the utmost importance, in order to 
ascertain possible thermal or engineering origins of the anomalous acceleration. 

The possibility that the gravitational pull from the Kuiper Belt may originate the 
anomalous acceleration was investigated in Refs. [1,7,8]; it was found that the reported 
magnitude of the effect would require that this extended object has a mass two orders 
of magnitude higher than the commonly considered value of Mxuiper — 0.3ME a rth- 

The inability to account for the Pioneer anomaly within the realm of conventional 
physics (see Ref. [9-11] and references therein) has driven several authors to consider 
a number of new theoretical ideas to explain it. However, proposals arising from new 
physics must be taken with due care, since an unambiguous description of the anomaly 
is still lacking. Indeed, the originally available Doppler measurements were conducted 
at distances that do not allow for a clear distinction of the direction of the acceleration: 
in particular, it is not known if the anomalous acceleration points towards the Sun or 
the Earth, along the line of sight. This information is of paramount importance: an 
effect pointing towards the Sun would indicate a gravitational origin (having discarded 
the far too weak hypothesis of solar wind pressure), while an anomaly directed at our 
planet would indicate that its origin lies in a modified Doppler effect (perhaps reflecting 
unknown physics affecting light propagation and causing an effective blue shift), or 
the incorrect interpretation of Doppler data, possibly due to systematic errors such 
as mismodeled Earth orientation parameters, incorrect ephemerides estimates, Deep 
Space Network or software clock drifts. Also, an anomalous acceleration along the 
spin axis of the spacecraft would show that it is caused by underestimated onboard 
systematic effects. 

If the acceleration is parallel to the velocity vector, this would indicate that the 
anomaly is a drag effect - although this does not appear to be physically motivated, 
as the reported magnitude of the acceleration requires that the probes fly through a 
medium with a density of order 10~ 19 g/cm 3 (see, e.g., Ref. [7]) - which should be 
compared with the density of interplanetary dust (arising from hot- wind plasma [13]) 
< 10 -24 g/cm 3 , and interstellar dust (directly measured by the Ulysses spacecraft) 
10~ 26 g/cm 3 . From a theoretical perspective, a modification of geodetical motion aris- 
ing from extensions of General Relativity could also account for a velocity dependent 
anomalous acceleration [14, 15] . 
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Initially regarded with some suspicion, the interest in the Pioneer anomaly has 
steadily grown, as disclosed by the number of peer-reviewed publications. The char- 
acterization of this anomalous acceleration was a major component of the scientific 
objectives of two mission proposals put forward to the recent ESA Cosmic Vision 
2015-2025 program [16]; although not considered by the reviewing board, this remains 
a burning question that may hint at yet unknown aspects of gravitational physics or, 
at the very least, allow us to increase our understanding of spaceflight dynamics. 



1.2 Previous Work 



Table 1 Error budget for the Pioneer 10 and 11, taken from Ref. [1]. 



Item 


Description of error budget constituents 


Bias 

10 -8 cm/s 2 


Uncertainty 
10" 8 cm/s 2 


1 


Systematics generated externa] to the spacecraft: 








a) Solar radiation pressure and mass 


+0.03 


±0.01 




b) Solar wind 




± < 10~ 5 




c) Solar corona 




±0.02 




d) Electro-magnetic Lorentz forces 




± < 10~ 4 




c) Influence of the Kuipcr belt's gravity 




±0.03 




f) Influence of the Earth orientation 




±0.001 




g) Mechanical and phase stability of DSN antennae 




± < 0.001 




h) Phase stability and clocks 




± < 0.001 




i) DSN station location 




± < 10~ 5 




j) Troposphere and ionosphere 




± < 0.001 


2 


On-board generated systematics: 








a) Radio beam reaction force 


+1.10 


±0.11 




b) RTG heat reflected off the craft 


-0.55 


±0.55 




c) Differential cmissivity of the RTGs 




±0.85 




d) Non-isotropic radiative cooling of the spacecraft 




±0.48 




c) Expelled Helium produced within the RTGs 


+0.15 


±0.16 




f) Gas leakage 




±0.56 




g) Variation between spacecraft determinations 


+0.17 


±0.17 


3 


Computational systematics: 

a) Numerical stability of least-squares estimation 

b) Accuracy of consistency/model tests 

c) Mismodclling of maneuvers 

d) Mismodclling of the solar corona 

e) Annual/diurnal terms 




±0.02 
±0.13 
±0.01 
±0.02 
±0.32 




Estimate of total bias/error 


+0.90 


±1.33 



Table 1, extracted from Ref. [1], provides an assessment of the systematic contri- 
butions to the acceleration budget; the different orders of magnitude of the various 
effects involved indicate that these do not account for the reported anomaly; how- 
ever, in another estimate of the heat dissipation of several spacecraft components, it is 
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claimed that a combination of sources could explain the anomalous acceleration [3]. A 
more recent and complete work has tackled the complex issue of modelling the Pioneer 
probes with great detail, aiming at the description of all relevant thermal effects with 
a sufficient accuracy [5]; parallel efforts are being pursued by other groups within the 
Pioneer collaboration team [6]. 

Preliminary results appear to indicate that up to one third of the total magnitude 
of the reported anomaly may be due to thermal effects [17]. However, it should be 
taken into account that the several modelling strategies and parameter estimations 
might cloud the overall picture - so that the physical significance could be difficult 
to distinguish amongst the technical depth of the computational reconstitution. For 
this reason, in Ref. [6] a complementary methodology is developed, that purposely 
diverts its attention away from the full modelling of intricate engineering detail, but 
instead focuses on the physical basis of the discussed thermal behaviour [6]. This ap- 
proach specifically assumes that, although a simplified modelling of specific details 
might diminish the global confidence of the obtained results, the added simplicity and 
computational promptness and speed is highly advantageous for a sensitivity analysis 
of the several relevant parameters. 

In this paper, the main features and results obtained so far with this method based 
on point-like sources are reviewed, including discussion of test cases, inclusion of diffu- 
sive reflectivity (and, in the future, specular), a preliminary order of magnitude budget 
for the thermal contributions from the spacecraft components, and its compatibility 
with previous studies. 



2 Source Distribution Method 

2.1 Motivation 

As discussed above, a thorough characterization of the Pioneer anomaly is required 
before more definitive statements about its origin can be made. This is the rationale 
behind the task of recovering and analyzing the full flight data, as well as the under- 
standing of the overall thermal behaviour of the Pioneer probes, which can give rise 
to a non-negligible acceleration. The key physical problem is to ascertain how thermal 
radiation is emitted and reflected by the external surfaces of the probes, and then 
calculate the resulting reaction force. 

Hence, alongside with the proposed finite elements model, a faster, more versa- 
tile approach was put forward [6]. This methodology, based on a distribution of a few 
point-like thermal sources, simulates the thermal radiation emitted from the spacecraft 
and takes into account what fraction is radiated directly into space, or reflected and 
absorbed by another surface of the spacecraft. This approach allows for a faster dis- 
crimination of different contributions from the major constituents of the vehicles: as a 
first approach, the RTGs, antenna dish and main bus compartment. 

In favour of this approach, one can state that it is impossible to model the Pioneer 
spacecraft in a very precise way: as it was built decades ago, the accuracy of the 
blueprints and existing models is limited; furthermore, such a long period in space 
leaves us with no detailed knowledge of the precise material properties, its degradation 
or damage. Hence, educated guesses will have to be taken, setting an intrinsic limit to 
the accuracy of obtainable results. 
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Furthermore, the fitting of thermal models to temperature data is also restricted 
by the data set (stemming from the onboard temperature sensors — six on the main 
bus and two on the RTGs) and lack of knowledge regarding the optical properties of 
the materials also yield uncertainties in the final result. Hence, it is clear that the total 
electrical power — which is well known — must be the fundamental parameter for any 
analysis. Our approach is based on this principle, besides its computational simplicity 
and focus on the relevant physics. 

Clearly, the distribution of power over the external surfaces of the spacecraft is 
essential to the quantification of thermal radiation effects. The insulation of the space- 
craft walls should limit the related gradient of the temperature along the main exter- 
nal surfaces (except in particular places — e.g. the louvers, which can be modeled as 
separate sources): we argue that small details and gradients in temperature of these 
external surfaces do not affect the results to a large extent, as shall be demonstrated 
by an analysis of the overall effect obtained upon variation of the number of point-like 
sources (keeping the power constant) . After considering these assumptions, a sensitiv- 
ity analysis of the several assumptions may be performed — including shape modelling, 
temperature gradients, and total power emitted by each surface; this sensitivity anal- 
ysis is advantageously achieved for the proposed model, due to its short computation 
time and straightforward physical basis. 

An important simplification relies on the spin stabilization of the Pioneer space- 
crafts, as the off-axis reaction forces will cancel out over time. Due to the geometry 
of the probes, the majority of small contributions possibly not taken into account are 
expected to be normal to the axis of rotation, and therefore irrelevant. This can be 
verified through a variation of the radiation distributions, which lead to similar values 
of the anomalous acceleration, as expected (however, the total reaction torque due to 
thermal radiation cannot be reliably predicted by the model here discussed). 

As hinted above, the simplicity of the formulation allows for a enhanced visibility of 
the involved physics throughout the entire modelling, thus allowing for a closer scrutiny 
of the method and derived results. Again, it is emphasized that the primary goal of our 
approach is to study a wide spectrum of the parameter space relevant for the thermal 
modelling of the Pioneer probes. This method is clearly not as comprehensive as a 
finite element model; however, it allows for a direct interpretation of results, higher 
adaptability and drastically reduced computation times. 

Before the key issue of the hypothetical thermal origin of the Pioneer anomaly is 
addressed, the self-consistency of the method should be strongly established: this is 
achieved through a battery of test cases that rely on simplistic geometrical arrange- 
ments. The validity of the adopted point-like source approach is also verified by a study 
of the convergence of results with the increase of the number of sources. 



2.2 Physical Formulation 

The method outlined in Ref. [6] is based on a distribution of isotropic and point-like 
sources so that, if W is the emitted power, the time-averaged Poynting vector for an 
isotropic source positioned at (xq, Ho, z o) is given by 

s . = YL (x - xp,y -yg,z - z a ) 

4yr [{x - xo) 2 + (y - yo) 2 + (z - z ff 2 
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A Lambertian source is characterized by a radiation intensity proportional to the cosine 
of the angle with the normal unit vector, 

s L b = _ cos 6 ( x - x o,y -yo, z - z o) / 2 n 

[(x - x ) 2 + (y - yo) 2 + {z - zo) 2 f /2 
Typically, isotropic sources are used to model point-like emitters and Lambertian 
sources to model surfaces. The resulting force and the amount of energy illuminat- 
ing a surface are obtained by integrating the Poynting vector of the source distribution 
over this surface. The total energy is given by the time-averaged Poynting vector flux, 

E ilum = Js-ndA = Js(G(s,t))- x dsdt . (3) 

where G(s, t) parameterizes the surface. Integration of the resulting force - as given by 
the radiation pressure multiplied by the unitary normal vector, yields the total force 
acting upon that surface. The radiation pressure is given by 

Prad = -S ■ n , (4) 

c 

where a radiation pressure coefficient 1 < a < 2 is introduced, so that a = 1 indicates 
full absorption, while a < 2 gives full diffusive reflection. 

The force acting on the source of the radiation may be calculated by integrating 
the radiation pressure multiplied by a normalized radial vector along a generic surface, 



/ 



Ml iA ■ (5) 

In an object endowed with a complex geometry (such as the Pioneer spacecrafts) 
shadows cast by the surfaces that absorb and reflect the radiation must be taken into 
account: the magnitude of this shadowing effect is obtained through a similar expression 
and then subtracted to the force obtained for the emitting surface. Alternatively, one 
may use an integration surface that includes the illuminated surfaces: the total result 
is the sum of all effects, a tota j = Fj/m. 



2.3 Test cases 

In order to demonstrate the efficiency of the method proposed in Ref. [6], a set of test 
cases is examined. The main question answered through this assessment concerns the 
ability to adequately represent thermal radiation emitted from an extended surface as 
an equivalent small number of point-like sources (opposed to considering many small 
thermal radiating elements). The test cases consider a square emitting surface with 
1 m 2 , so that the force on the emitting surface, shadow caused by another surface at 
a given position and radiation pressure on that surface may be computed. The results 
for different numbers of sources are compared, keeping the total power constant: the 
result should converge to the exact solution with the increase of the number of radiation 
sources. In Ref. [6] it is shown that a reliable error estimate may be obtained modelling 
the surface with just a small number of sources. 

Taking a distribution of Lambertian sources on a surface on the Oxy plane and 
using Eq. (5) one obtains, for radiation emitting surfaces without any other illuminated 
surfaces, a force pointing in the z-axis and of magnitude (2/3)W sur f /c, depending only 
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I (W sr 1 ) 




Fig. 1 Polar plot of the intensity variation with elevation of the radiation emitted by a surface, 
when considering 1, 4, 16, 64 or 144 Lambcrtian sources, maintaining the total emitted power 
constant. The intensity diminishes with the number of sources. 



on the total emitted power. However, the calculations of the shadow and pressure 
radiation on other surfaces do depend on the source distribution: the variation of the 
radiation intensity with the elevation and the azimuth for 1, 4, 16, 64 and 144 source 
meshes is displayed in Figs. 1 and 2. These plots show that the maximum deviation 
occurs at the higher elevations and is less than 10%, for just one point-like source. 
Actual deviations for angles relevant to the Pioneer will be considerably smaller. This 
estimate may be confirmed by calculating the force acting on another 1 m 2 surface, for 
several different positions: a total of nine configurations were taken, each with different 
positions and tilt angles (see Table 2); the deviation between the 1, 4, 16, 64 and 144 
source meshes is then obtained. 

The highest deviation was obtained for test case 8, confirming the assumed scenario 
— since the second surface is at a high elevation from the emitting surface, as shown 
in Fig. 3. The results in Table 3 indicate a difference of about 6% between the force 
obtained with one source and finer meshes (16, 64 and 144 sources). The latter are 
within 0.5% of each other, and the intermediate 4 source mesh has a deviation of just 
1.5%. 

Test cases 1 and 3 are the most relevant, since these depict configurations typical of 
the Pioneer spacecraft: the first case yields results present in Table 4 (see also Fig. 4), 
showing that, for 16, 64 and 144 sources, variation in the energy flux and force is 
smaller than 0.5%. Furthermore, comparison with finer meshes shows that results vary 
by less than 5% for one source, and less than 1.5% for a four source mesh. For test case 
3, the results in Table 5 show a variation of less than 5% between the results for one 



s 



I(W rad" 1 ) 




Fig. 2 Same as Fig. 1 but for intensity variation with azimuth. Intensity diminishes with the 
number of sources. Notice that beyond 16 sources the results are almost coincident. 

y 




Fig. 3 Geometry of test case 8 (cf. Table 2): thermal emission from a surface is simulated by 
a varying number of Lambertian sources evenly distributed on the surface, keeping the total 
emitted power fixed, and the effect on another surface is observed. This test case leads to the 
highest variation of results with the number os sources. 
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Table 2 Positions considered for the second surface in test cases; distances between both 
surfaces are typical for the Pioneer spacecraft. 



Test case 
# 


Surface Centre Position 
(m) 


Surface Tilt Angle 
(°) 


1 


(0, 2, 0.5) 


90 


2 


(2, 0, 1.5) 





3 


(2, 0, 1.5) 


.30 


1 


(2, 0, 1.5) 


60 


5 


(2, 0, 1.5) 


90 


6 


(1, 0, 2) 





7 


(1, 0, 2) 


.30 


8 


(1, 0, 2) 


60 


9 


(1, 0, 2) 


90 



Table 3 Results for test case 8 (cf. Table 2) considering a total emission power of 1 kW. The 
force components due to shadowing on the secondary surface are almost independent of the 
number of sources on the emitting surface. 



Sources 


Energy Flux 


Force components (x, y, z) 


# 


(W) 


(10- 7 N) 


1 


45.53 


(2.016, 0, 2.083) 


4 


45.53 


(1.918, 0, 2.003) 


16 


45.53 


(1.895, 0, 1.984) 


64 


45.53 


(1.890, 0, 1.979) 


144 


45.53 


(1.889, 0, 1.978) 




Fig. 4 Same as Fig. 3 for test case 1. 



and 144 sources. As in the previous cases, convergence is achieved for the 16, 64 and 
144 source meshes, with a variation of less than 0.25%. 

Ultimately, one concludes that a mesh with four Lambertian sources, which shows 
deviations of about 1.5%, provides the desired balance between precision and simplicity; 
this yields a convenient illustration of the ability of the developed method to estimate 
the radiation effects on the Pioneer probes. Given that the deviations are always well 
below 10%, even when the roughest approximations are considered, it appears that the 
application of this source distribution method to the scales and geometry involved in 
the Pioneer anomaly problem yields consistent and convergent estimates of the thermal 
radiation effects. 
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Table 4 Same as Table 3 for test case 1. 



Sources 


Energy Flux 


Force components (x, y, z) 


# 


(W) 


(10- 7 N) 


1 


15.34 


(0.9300, 0, 0.1514) 


4 


15.92 


(1.028, 0, 0.1638) 


16 


16.09 


(1.038, 0, 0.1675) 


64 


16.13 


(1.040, 0, 0.1684) 


144 


16.14 


(1.041, 0, 0.1686) 



Table 5 Same as Table 3 for test case 2. 



Sources 


Energy Flux 


Force components (x, y, z) 


# 


(W) 


(10- 7 N) 


1 


19.20 


(0.4952, 0, 1.037) 


4 


19.83 


(0.5032, 0, 1.082) 


16 


19.99 


(0.5050, 0, 1.093) 


64 


20.03 


(0.5054, 0, 1.096) 


144 


20.04 


(0.5055, 0, 1.096) 



3 Thermal Radiation Model of the Pioneer Spacecraft 



3.1 Geometry 

A simplified model of the Pioneer spacecraft can be obtained by considering some 
plausible, physically motivated hypotheses, such as the aforementioned probes' spin 
stabilization, and also the assumption of steady-state thermal equilibrium throughout 
most of their journey. Also, a simplified version of the spacecraft geometry may be 
considered as a first approach to the problem of estimating its thermal radiation: 
hence, one considers only its most important features (as illustrated in Fig. 5): the 
parabolic antenna, RTGs and a prismatic equipment compartment, with dimensions 
derived from the available Pioneer technical drawings. The estimation of thermal effects 
proceeds by integrating the emissions of the RTG and lateral walls of the equipment 
compartment along the visible portion of the antenna; since the emissions stemming 
from the front surface of the Pioneer probes is never reflected, this element is considered 
as a whole. Furthermore, one may discard the surface of the compartment that faces the 
antenna, since its contribution to the anomalous acceleration is negligible for obvious 
geometric reasons: its emitted radiation, which is mainly radial, is attenuated through 
multiple reflections. Finally, the antenna should have a very low, approximately uniform 
temperature (as discussed in Ref. [17]), also yielding a negligible effect. 

It will be demonstrated how this simplified model is able to reproduce the most 
important contributions to the thermal reaction force: the RTGs and the main equip- 
ment compartment are effectively responsible for the overwhelming part of the emitted 
thermal radiation. In the latter case, most of the effect arises from the louvers located 
in the front wall (facing away from the Sun) - with consequences for the total power 
distribution. 
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Fig. 5 Pioneer spacecraft model geometry, back view: high gain parabolic antenna and hexag- 
onal main bus compartment. 

3.2 Point-like Source Distribution 

An analysis of each of the three main contributions to the total thermal effects is 
now developed. The front wall of the spacecraft, where an axial force with magnitude 
(2/3)Wf ront /c. The side walls of the main compartment is calculated by integrating the 
shadow and radiation pressure components along the antenna. One may neglect the 
shadow of the RTGs, due to their smallness and relatively distant positioning, which 
would yield a contribution mostly in the radial direction. 

Following the approach depicted in the test cases already discussed, the integration 
is performed with an increasing number of sources, thus allowing for evaluation of the 
method's convergence: this occurs rather quickly, with deviations consistently below 
2.5%. The values obtained indicate that 16.8% to 17.3% of the power emitted from the 
side walls of this compartment is converted into approximately Sun-ward thrust along 
the z-axis. 

The relevance of assuming a non-uniform temperature distribution is also assessed, 
by varying the relative power of the point-like sources in each surface, keeping the total 
surface power constant. A variation of 20% in power between sources, which models a 
5% temperature variation, does not produce significant changes in the final result — 
since observed relative differences are smaller than 1%. 

Finally, the RTG contribution is computed through two different models. The first, 
simpler scenario, mimics each RTG with a single isotropic source. In this case, the 
point-like source has the whole power of the RTG. In the second model, the cylindric 
shape of the RTG is taken into account and a Lambertian source is placed at each 
base. Actually, it is only necessary to consider the source facing towards the centre 
of the spacecraft, as the remaining RTG radiation will be emitted radially and its 
time-averaged contribution vanishes. In this case, the Lambertian source has a certain 
amount of the total RTG power, as discussed in the following Sections. Depending on 
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the model considered, either 1.9% of the total power or 12.7% of the power emitted from 
the base of the cylinder (equivalent to approx. 2% of total RTG power, if temperature 
is uniform) is converted into thrust. 

In the following section, more accurate figures will be obtained, which include the 
diffusive reflection included in Eq. (4) . 



3.3 Available Power 

Amongst the panoply of physically meaningful parameters, the available power on the 
Pioneer spacecraft is one of the few that is reasonably well known. For this reason, 
the methodology pursued here takes this quantity as the independent variable from 
which all estimates of the resulting thermal effects are derived, instead of temperature 
readings - which must be interpolated from the isolated sensor sites to the remainder 
elements of the spacecraft, an objective outside of the scope of this study. Clearly, the 
(assumed uniform) temperature of a surface i is related to the energy balance to the 
spacecraft in steady-state conditions, through Stefan's law, 



where Ai are the relevant areas and ej the emissivity of surface i. 

Since one does not have accurate information about the optical properties of the 
surfaces, as well as their evolution in time, the temperature estimates are quite uncer- 
tain. Also, note that a variation in the emissivity requires the computation of a new 
solution for the temperature distribution, in order to obtain the correct total power 
and maintain energy conservation. 

The power generated onboard the probes comes from the two Plutonium-238 pow- 
ered RTGs, with just a fraction of the generated heat converted into electricity: the 
remaining power is dissipated as heat, which will be transfered to the central com- 
partment through conduction in the truss assembly. However, given the small section 
of this structure, one may disregard this contribution and take the RTGs as isolated 
elements, with all produced thermal power radiated directly from them. 

A considerable proportion of the electrical power is consumed by the high gain 
antenna, with the remainder being used by the various instruments located in the main 
compartment. As found in Ref. [1], the total RTG thermal power at launch amounts to 
2580 W, yielding 160 W of electrical power - so that, at launch, approximately 2420 W 
of thermal power was radiated away by the RTGs. Considering the plutonium decay 
with a half-life of 87.74 years, the total on-board power variation with time (in years) 
is given by 



However, telemetry data reveals that the electrical power decayed at a faster rate: the 
electrical heat in the body of the spacecraft was around 120 W at launch, dropping 
to less than 60 W at the latest stages of the mission [5], following an approximate 
exponential decay with a half-life of only 24 years. This accelerated power decay reflects 
the added effect of thermocouple degradation, which reduces the conversion efficiency 
of the RTGs. 




(6) 




(7) 
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4 Results and Discussion 

4.1 Order of Magnitude Analysis 

Before proceeding to a more rigorous numerical estimate, the results detailed in the 
previous Sections can be used to perform a preliminary order of magnitude analysis; 
this will serve to produce a concrete figure of merit for the overall acceleration arising 
from thermal effects, which can then be compared with the dp ~ 10 -9 m/s 2 scale of 
the Pioneer anomaly. 

The spacecraft specifications indicate a total mass for the probes mpi ~ 230 kg, 
and RTG and equipment compartment powers Wrtg ~ 2 kW and VK e quip ~ 100 W, 
respectively. As discussed, integrating the corresponding thermal emissions yields the 
proportion of emitted power that is effectively converted into thrust. Assuming a uni- 
form temperature and emissivity in the RTGs and main equipment compartment, which 
leads to a power emitted from each surface proportional to its area; with the simple 
model discussed is Section 3.2, this yields 

F RTG ~2xl0- 2 ^^, (8) 

-1 W^equip -1 Wcquip 

^sidcs ~ 10 " , front ~ 2 x 1U " ■ 

c c 

The acceleration of the spacecraft due to the thermal effects arising from the power 
dissipation of the RTGs and equipment compartment is easily obtained, 



CJRTG ~2X 10 -2liRTG _ 6 x 10 -10 m/g 2 (g) 
m Pio c 



mpi c 



1 c w in— 1 "equip no., in— 10 / 2 

floquip ~ 1-5 x 10 m ~ 2.2 x 10 m/s 



so that both contributions are relevant when estimating the acceleration induced by 
thermal radiation from the Pioneer probes: since the RTGs and the instrument com- 
partment produce similar thermal effects, one cannot rely solely on one of these sources, 
confirming the previous discussion in Refs. [3,5]). 



4.2 Thermal Force Estimate 

Having gained momentum through the estimate depicted above, a more thorough eval- 
uation of the existing thermal effects using the point-like source modelling developed 
in Ref. [6] follows: to this effect, one considers a model with 4 sources per side panel 
of the equipment compartment, and Lambertian sources at the bases of the RTGs (see 
discussion in section 3.2) - the configuration that provides what is deemed as the best 
compromise between accuracy and computation time, with the aforementioned devi- 
ation < 0.5% from results using finer meshes. The acceleration is obtained from the 
axial component of the integration of the radiative emissions of the Pioneer probes, 

(0.168H/ Sidcs + §VK front + 0.128M/ baso ) 

apio = > (10) 

m Pio c 
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Fig. 6 Variation of the resulting acceleration with the temperature ratio between the louvers 
and the equipment platform, considering similar cmissivities for both multi-layer insulations. 



where W s id C s an d W{ roni are the powers emitted from the side panels and front of the 
equipment compartment and IVbasc is the power emitted from the base of the RTG 
facing the centre of the spacecraft. 

To ascribe values to each of these powers, one takes readings from 1998, as found in 
Ref . [5] , namely the dissipated thermal powers at the RTG and equipment compartment 
Wrtg = 2050 W and W equ i p = 58 W, respectively. With the simplest scenario of 
uniform temperature and optical properties, this gives 



Wsides = 21.75 W , VK front = 18.12 W (11) 
W hase = 41.11 W , 

leading to an acceleration a^h = 3.05 x 10~ 10 m/s 2 . This amounts to about 35 % of 
the anomalous acceleration. 

From the available temperature maps used in Refs. [5,17], one can see that the 
temperature anisotropies along the sides of the equipment compartment rest within 
the limits assumed in the test cases, as discussed in Section 3.2. However, the RTGs 
display significant temperature variation between the wall of the cylinder, the bases 
and the fins. Also, the front wall of the equipment compartment is expected to produce 
a larger contribution than that of the side walls, due to the presence of the louvers. 

This said, one can ascertain how does the emitted power in the louvers and at 
the base of the RTG vary — two critical parameters in this calculation. If the louvers 
are closed and assumed to have a emissivity similar to the equipment platform, the 
variation of the acceleration with the temperature ratio between the louvers and the 
mean temperature of the platform may be plotted (again, with the total power fixed) , 
as shown in Fig. 6. The same analysis can be performed with respect to the RTGs, 
now considering the ratio between the temperatures at the base of the cylinder and 
the fins (shown in Fig. 7) . 

Figs. 6 and 7 help to better grasp the main advantage of the developed method 
[6]: the possibility of conducting a sensitivity analysis that is both fairly quick and 
sufficiently accurate; from Eq. (10) and a plausible variation of the power parame- 
ters, temperature readings may be adjusted and different optical properties suitably 
discussed and inferred. 
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Fig. 7 Variation of the resulting acceleration with the temperature ratio between the base of 
the RTG cylinder and the fin temperature. 

Another scenario that may be explored assumes that the RTG cylinder bases and 
wall possess a temperature 15% and 30% higher than that of the fins, respectively, 
while the closed louvers have double the temperature of the equipment compartment 
(with constant emissivity). This yields the following estimates for the power, 

W sidcs = 9.97 W , WW = 39.71 W , (12) 
WW = 49.67 W , 

thus accounting for 57% of the anomalous acceleration, that is, aj~h = 5.00x 10~ 10 m/s . 

Up to this point, the study has considered only the case of full absorption of the 
radiation by the illuminated surfaces; to counteract this, one includes diffusive reflec- 
tion by assigning a non-unitary value to the a parameter in Eq. (4). For the relevant 
wavelengths, the reflectivity of the type of aluminum used in the high gain antenna is 
commonly around 80% , yielding a = 1.8; also, the multi-layer insulation of the equip- 
ment platform may be modeled by a value a — 1.7. Given these conditions and the 
re-derived illumination factors in Eq. (10), the same temperature conditions considered 
in the previous case result in a thermally induced acceleration a Th = 5.75 x 10~ 10 m/s 2 . 

Besides allowing for the computation of observables, the derivations put forward in 
this Section assess the variations occurring when considering different parameters and 
hypotheses. The three discussed scenarios serve to show that the method developed 
in Ref. [6] is well suited to identify sensitive parameters and promptly evaluate the 
propagation of the existing uncertainties. 

5 Conclusions 

In this contribution we discuss the applicability of the method developed in Ref. [6] 
to account for the acceleration of the Pioneer spacecrafts arising from thermal effects, 
based on point-like sources. After identifying the main power contributions of the var- 
ious components of the spacecraft, 35% to 57% of the reported anomalous acceleration 
is obtained. 
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The discussed method displays a reasonable degree of accuracy, with a direct nu- 
merical computation error of the order of 1CP 1 or less, and deviations due to the 
approximation of the geometry with point-like sources falling below 1% (as argued in 
Sections 2.3 and 3.2). However, it should be emphasized that this is not a direct in- 
dication of the accuracy of the actual thermally induced acceleration, when compared 
to the reported case of the Pioneer anomaly — but a measure of self-consistency of 
the developed method. The aim here is to demonstrate the reliability of the proce- 
dure, which should be extended in order to better simulate the physical system of the 
Pioneer spacecrafts - while keeping the flexibility and computational speed that have 
motivated this complementary approach. 

Thus, future refinements will focus on an enhanced geometrical modelling, including 
specular reflection. Furthermore, it will be of great interest to identify the parameters 
that most directly affect the final result, e.g. temperatures, emissivities and reflectivities 
of the various components. Clearly, the encouraging results of the analysis performed 
in Ref. [6] are intimately connected to the level of agreement with thermal models 
based on finite element methods: as this work demonstrated, our estimates are so far 
consistent with these approaches. 
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